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ABSTRACT 

We report an over-density of bright sub- millimetre galaxies (SMGs) in the 0.15 deg 2 
AzTEC/COSMOS survey and a spatial correlation between the SMGs and the optical- 
IR galaxy density at z Si 1.1. This portion of the COSMOS field shows a ~ 3tr over- 
density of robust SMG detections when compared to a background, or "blank-field" , 
population model that is consistent with SMG surveys of fields with no extragalactic 
bias. The SMG over-density is most significant in the number of very bright detections 
(14 sources with measured fluxes Si.i mm > 6 mJy), which is entirely incompatible with 
sample variance within our adopted blank-field number densities and infers an over- 
density significance of » 4cr. We find that the over-density and spatial correlation 
to optical-IR galaxy density are most consistent with lensing of a background SMG 
population by foreground mass structures along the line of sight, rather than physical 
association of the SMGs with the z < 1.1 galaxies/clusters. The SMG positions are 
only weakly correlated with weak- lensing maps , suggesting that the dominant sources 
of correlation are individual galaxies and the more tenuous structures in the region 
and not the massive and compact clusters. These results highlight the important roles 
cosmic variance and large-scale structure can play in the study of SMGs. 

Key words: surveys - galaxies: evolution - cosmology: miscellaneous - infrared: 
galaxies submillimctcr gravitational lensing 



1 INTRODUCTION 

Foreground structure, cosmic variance, and source environ- 
ment can affect the observer's perception and interpreta- 
tion of the source population being probed in a particular 
survey. For example, gravitational lensing by massive fore- 
ground clusters affects both the observed flux of sources 
and the areal coverage of the survey in the source plane. 
These aspects of gravitational lensing have been utilised to 
probe the very faint sub-millimetre galaxy (SMG) popula- 
tion below the confusion limit imposed by the high den- 
sity of faint SMGs relative to the survey beam size (e.g . 
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Il997l ; IChapman et al.l l200d; ICowie et al 
[2002! : iKnudsen et al.ll2006l ; fwilson et al ' 



2002; 



2008 



The measured (sub) millimetre fluxes of sources found in the 
direction of very massive clusters can also be affected by, 



and confused with, the signal imposed through the Sunyaev- 
Zel'dovich effect o n the cosmic microwave background (e.g. 
I Wilson et al.|[2008h . Furthermore, surveys can be affected by 
foreground structures with high galaxy-densities, which in- 
crease the likelihood of galaxy-galaxy lensing by intervening 
galaxies and complicate counterpa r t identification at o ther 
wavelengths (IChapman et al.ll2002l : iDunlop et al.ll2004 ). 

Spectroscopic observations have shown that the vast 
majority of SMGs with detect able radio counterpart s lie at 
an average red shift of z ~ 2.2 (IChapman et al]|2005l ). while 
spectroscopic ijValiante et all 120071 ) and photometric (e.g. 
lYounger et a.1 .1120071 ') analysis put many radio-faint SMGs at 
even higher redshifts. The average SMG is unlikely to be 
found at z <, 1, however it remains to be seen if the z ~ 1 
SMG population can be locally enhanced due to large-scale 
structure and cosmic variance. Some evidence exists for in- 
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creased number densities of SMGs in mass-biased regions 
of th e z > 1 Universe. Surveys to wards several z ~ 1 clus- 
ters (|Bestll2002l : IWebb et al]l2005h find a number density of 
SMGs in excess of the blank-field counts that can not be 
explained by gravitational lensing alone. This implies that 
some of the SMGs are physically associated with the clus- 
ters, although the number statistics are small and the lensing 
could be underestimated l|Webb et al.l l2005h . Similar over- 
den sities have been found towards high-redshift radio galax- 
ies (IStevens et al.ll2003l : |Pe Breuck et al J 12004 iGreve et all 
120071 ') and z > 5 quasars (jPriddev et al where lensing 

of background sources is less likely to be an issue. Spec- 
troscopic observations have also found common redshifts 
amongst SMGs in the SSA22 and HDF fields, suggesting 
physical over -densities of SMGs at redshifts of 3.1 and 2.0, 
respectively (|Chapman et alj|2005h . Together, these surveys 
suggest that these massive dusty starbursts are prominent 
in moderate and high-redshift cluster/proto-cluster environ- 
ments. 

In this paper, we analyse the density a nd distribu- 
tion of SMGs in the AzTEC/COSMOS survey (IScott et all 
120081 ). The AzTEC/C OSMOS survey cove rs a region within 
the COSMOS field l|Scoville et all 120071 ) known to con- 
tain a high density of optical-IR galaxies and prom inent 
large-scale structure at z <> 1.1 (|Scoville et al.l I2007T ). in- 
cluding a massive M ~ 10 15 M Q cluster at z « 0.73 
l|Guzzo et al.l I2007T ). In § [5] we present the 1.1 mm source 
counts for the AzTEC/COSMOS field, revealing a strong 
over-density of bright SMGs compared to the blank-field. 
We explore the nature of this over-density through exam- 
ination of the spatial correlation between SMGs and the 
known large-scale structures (§[3|. Positive correlation be- 
tween SMG positions and low-redshift large-scale struc- 
ture has been previo usly detected statistically in 3 dis- 
joint SCUBA surveys (|Almaini et al.ll2003l . h00$ ). We now 
present a wide-field investigation of such correlations using 
the AzTEC/COSMOS survey, which has advantages in its 
contiguous size (0.15 deg 2 ), broad range of low-redshift en- 
vironments, and the availability of deep multiband im aging 
and reliable photometric redshifts jllbert et al.ll2008r ). 

This is the second paper describing t he 1.1 mm results 
of the AzTEC/COSMOS survey. Paper I (|Scott et al.ll2008l ) 
presented the data-reduction algorithms, AzTEC/COSMOS 
map and source catalogue, and confirmation of robust- 
ness of the AzTEC/JCMT data and pointing. Addition- 
ally, seven of the brightest AzTEC/COSMOS sources have 
had high-resolution follow-up imaging at 890 [im using the 
Sub-Millimetre Arra y (SMA) and are discussed in detail in 
lYounger etHI (|2007l ). Spitz er IRAC co l ours o f these SMGs, 
and others, are discussed in lYun et alj l|2008f ). 



2 NUMBER COUNTS IN THE COSMOS FIELD 

The number density of SMGs provides const r aints on galaxy 
evolution models (e.g. | Kaviani et al.l 20031; iGranato et al.l 



12004 iBaugh et all 120051 : iNegrello et all 120071 ') and insights 
to the dust-obscured component of star formation in the 
high-redshift Universe. The number density also describes 
how these discrete objects contribute to the cosmic infrared 
background (CIB), as discussed in Paper I. In this paper, we 
focus on how the localised SMG number counts reflect large- 



scale structure. Before presenting the number counts for the 
AzTEC/COSMOS survey (§H3]& §[23]), we describe the 
technical details of the flux corrections (§ I2.1|l and methods 
(§ 12. 2[) that are vital to the construction of unbiased source 
counts from typical SMG su rveys. Here w e expa nd on the 
flux correction techniques of ICoppin et all (|2006l ) and pro- 
vide new tests of these methods through simulation. 



2.1 Flux Corrections 

Surveys of source populations whose numbers decline with 
increasing flux result in blind detections that are biased sys- 
tematically high in flux. This bias is typically referred to 
as "flux boosting" and results from the fact that detected 
sources have a higher probability of being an intrinsically 
dim source (numerous) coincident with a positive noise fluc- 
tuation than being a relatively bright source (scarce) coin- 
ci dent with negative noise . This effect is concisely described 
in iHogg &: Turned (|l998T ) and is extremely important for 
SMG surveys (see Figure [1} due to the relatively low S/N 
of the measurements and a population that is known to 
decline steeply wit h increasing flux (e.g. IScott et al. I l2006l : 
ICoppin et al.l 20061 . and references therein) . 

We calculate an intrinsic flux probability distribution 
for each potential AzTEC sour ce using the Bayesian tech - 
niques outlined in Paper I and ICoppin et al. (2005, 2006). 
The probability of a source having intrinsic flux Si when 
discovered in a blind survey with measured flux S m ± <r m is 
approximated as 



p(S t \S„ 



p(Si)p(S m ,a m \Si) 

p(S m , 0"m) 



(1) 



where p(Si) is the assumed prior distribution of flux densi- 
ties, p(S m , <Jm\Si) is the likelihood of observing (S m , o m ) for 
a source of intrinsic flux Si, and p(S m , <y m ) is a normalising 
constant. The resulting probability distribution is referred 
to as the posterior flux distribution, or PFD, throughout 
this section. We assume a Gaussian (normal) noise distri- 
bution for p(S m , a m \Si) that is consistent with the noise in 
our map at the location of the discovered source. The prior, 
p(Si), is generated from pixel histograms of 10,000 noiseless 
simulations of the astronomical sky - as would be seen in 
zero- mean AzTEC/JCMT maps - given our best estimate of 
the true underlying SMG population and distribution. For 
this paper, we assume the SMG population exhibits num- 
ber coun t densities that are well described by a Schechter 
function (Schechter 1976) of the form 



dN 
~dS 



AT 
~S- 



T (J) exp(-5/S'). 



(2) 



This parametric form is a slight departure fro m that used in 
Pape r I and in the SCUBA/SHADES survey (jCoppin et all 
2006), with N* / S' replacing the parameter N' found in Pa- 
per I. This form has the advantage of reducing the cor- 
relations between the normalising parameter and the pa- 
rameters S' and a. The normalising factor, N* , is in units 
of deg -2 and is independent of the observation wavelength 
when assuming the same source population and a constant 
flux ratio between the observing bands. For the Bayesian 
prior, we initially assume parameters of [S',N*,a] — 
[1.34, 5280, —2], which represent the best-fit S chechter func- 
tion to the SCUBA/SHADES number counts (jCoppin et all 
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2006) when reparameterised to the form of Equation [2] and 
scaled to 1.1 mm assuming an 850 u.m/1100 u.m spectral 
index of 3.5 (flux ratio ~ 2.5). 

A second systematic flux bias in low-S/N blind sur- 
veys results from source detections being defined as peak 
locations in the map. The measured source flux is, on av- 
erage, biased high due to the possibility of large positive 
noise peaks lying nearby, but off-centre from, the true source 
position. This bias is minimised through point-source fil- 
tering and is sub-dominant to the flux boosting described 
previously. It is significant only for the lowest S/N detec- 
tions and is largely avoided by restricting our analysis to 
the most robust sources (S/N 4). The remaining small 
bias (bpeak < 0.2a m for S/N ^ 4) is estimated through sim- 
ulation and subtracted from the detected source flux (S m ) 
before calculating the PFD in Equation [1] 

We validate these flux corrections through extensive 
simulation of the PFD. We generate 10,000 simulated maps 
by adding noiseless sky realisations to random noise maps 
using the prescription outlined in Paper I. Sources are ran- 
domly injected spatially (i.e. no clustering) and in accor- 
dance with the number counts prior assumed. We group 
recovered sources in the resulting maps according to their 
measured values (S m ,a m ), with each being mapped back 
to an intrinsic flux, Si, defined as the maximum input flux 
found within Obeam = 7.6" of the output source location. For 
each bin of measured values (S m ,crm), the input Si values 
are binned and normalised to produce a simulated PFD. 

Example simulation results are presented in Figure [1] 
Overall, the Bayesian approximation of the PFD (solid 
curve) provides a good estimate of the simulated proba- 
bility distribution (histogram) at most fluxes. The differ- 
ences at low flux and low S/N are due to a combination 
of source confusion in the simulations, higher-order effects 
of the bias to peak locations, and other low level system- 
atics. For the purposes of this paper, the Bayesian results 
are preferred to the simulated PFDs for their computa- 
tional speed, resolution, and flexibility in priors. The strong 
differences between the simulated probability distributions 
(histograms) and the naive Gaussian distributions (dashed 
curves) demonstrate the significance of flux boosting in sur- 
veys of this type. It is important to note that flux boost- 
ing (as described above) is not related to the adopted de- 
tection threshold and that even the most robust detections 
can be significantly biased. For example, a source detected 
at S/N = 8 in the AzTEC/COSMOS map will have been 
boosted by an average of 1.2 mjy (~ la m ), assuming the 
scaled SCUBA/SHADES SMG population. 

2.2 Number Counts Derivation 

The relative robustness of each source candidate is encoded 
in the PFD and is a function of both S m and a m , as op- 
posed to merely Sm/^m, due to the population's steep lu- 
minosity function. As in Paper I, we use the total prob- 
ability of a source candidate being de-booste d to negative 
flux as the metric of relative source robustness. ICoppin et al.l 
(|2006h found that P(S l < 0\S m , cr m ) < 0.05 provided a nat- 
ural threshold from which to select a large sample of ro- 
bust SMGs without including a significant number of noise 
peaks, or "false detections". This threshold also marks the 
point where the Bayesian approximation begins to suffer 
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Figure 1. Posterior Flux Distribution (PFD) for S/N = 4 (left) 
and S/N = 5 (right) detections in a map with noise ~ 1.3 mjy 
assuming an underlying source population consistent with the 
AzTEC/COSMOS results presented in §[2T3] The simulated prob- 
ability distribution is shown as the histogram and is calculated 
as described in § 12.11 Error bars represent the lcr Poisson errors 
of the simulation results, limited only by the number of simu- 
lations computed. The Bayesian approximation is depicted as a 
solid curve. The dashed vertical line represents the measured flux, 
S m - The dashed curve is the Gaussian probability distribution, 
p(Sm, o" m | Si), which represents the distribution that might oth- 
erwise be assumed without flux boosting and/or false detection 
considerations. 

from low-level systematics, as suggested by the comparison 
of the Bayesian and simulated PFDs (Figure [TJ. Therefore, 
we will use this "null threshold" of 5% to define our cat- 
alogue of robust sources from which to estimate number 
counts. This threshold is equivalent to S/N values of 4.1 - 
4.3 for our range of a m values, 1.2 - 1.4 mjy, assuming the 
scaled SCUBA/SHADES prior. 

We derive the number counts from the catalogue of 
robust sources and their associated PFDs using a boot- 
strap sampling method similar to that used in lCoppin et al.l 
(2006). In each step of this method, the selected sources 
are randomly assigned fluxes according to their respective 
PFDs (Equation [TJ). These samples are binned by flux to 
produce differential (dN/dS) and integral (N(> S)) source 
counts, with each bin being appropriately scaled for survey 
completeness and area. We introduce sample variance by 
sampling the rob ust source catalogue with replacement (e.g. 
iPress et al.lll992l ), and by Poisson deviating the number of 
times the catalogue is sampled around the true number of 
detections. We repeat this process 20,000 times to determine 
uncertainty and correlation estimates for the number count 
bins. 

Applying this sampling method to relatively small 
source catalogues results in a discretely sampled probability 
distribution for each number counts bin. This finite multi- 
nomial distribution can be non-Gaussian and asymmetric; 
therefore, we describe the uncertainty in the number counts 
as 68% confidence intervals that are approximated by lin- 
early interpolating between the occupation numbers sam- 
pled in the bootstrap. 

Survey completeness is estimated through simulation, in 
which sources of known intrinsic flux are randomly injected 
into noise map realisations one at a time and their output is 
tested against the null threshold source definition. Indepen- 
dent simulations confirm that this method provides excellent 
completeness estimations at all fluxes considered and that 
source confusion is not an issue given our beamsize (18") 
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and map depth (a > 1.1 mJy). Completeness is calculated 
as a function of intrinsic flux and averaged across the map 
to account for the slightly varying depth across the survey 
region considered. We calculate the effective completeness of 
each differential number counts bin by averaging the simu- 
lated completeness function within the bin, weighted by the 
assumed relative abundance of sources (i.e. the prior). 

We test these techniques by applying the same flux cor- 
rection and number counts extraction algorithms to sim- 
ulated maps with the same size and noise properties as 
those of the AzTEC/COSMOS survey. Figure [2] shows the 
extracted differential number counts from simulated maps 
using two different assumed priors. Both sets of simu- 
lated maps were populated with the same SMG popula- 
tion (solid line), which is similar to the final results of this 
AzTEC/COSMOS survey (§ E3J. The right panel of Fig- 
ure H] shows the results of the ideal case where the Bayesian 
prior is the same distribution used to randomly populate the 
simulated maps, while the left panel shows the results when 
using the scaled SCUBA/SHADES prior (dashed curve), 
which differs from the simulated input population (solid 
curve). For both priors the extracted number counts are in 
excellent agreement with the injected population. The rela- 
tively small differences between the input and output counts 
in the ideal case are used as systematic correction factors in 
our final calculations. The lowest flux bin (1-2 mJy) suffers 
from very low (and poorly defined) completeness and is, in 
general, the most sensitive to the assumptions in the prior. 
For these reasons we will restrict our analysis in this pa- 
per to the number count results for fluxes > 2 mJy, unless 
otherwise specified. 

In Figure [2] the dispersion of the output source counts 
(error bars) is noticeably smaller than the dispersion of in- 
put source counts (shaded region) at high fluxes. This dis- 
crepancy reflects the correlation between output data points 
through our assumed prior. We characterise the overall bias 
to the assumed population by testing a wide range of priors 
against a static input population. For priors that are consis- 
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Figure 2. Simulated differential number count results (data 
points) using the extraction techniques described in § I2.2I The 
solid line represents the analytic source counts (Equation [2j used 
to populate the simulated maps, while the shaded region is the 
1<t dispersion incurred by randomly populating maps of this size 
(0.15 deg 2 ). The left panel shows the results when assuming 
a prior based on the scaled SCUBA/SHADES results (dashed 
curve), while the right panel is for the "ideal" prior that matches 
the underlying input population (solid curve). Error bars repre- 
sent the ltr dispersion from 1,000 simulations, while the errors in 
the means are typically smaller than the data symbols plotted. 



tent with previous SMG surveys (e.g . lLaurent et al.l [20051 : 
ICoppin et ail 120061 ; | Scott et aHl2006t i. the bias incurred is 
generally smaller than the formal la errors of the extracted 
counts in a survey of this size and depth. Larger biases can 
result for exceptionally poor priors (e.g. > order of mag- 
nitude differences from the true population); however, in 
most cases the extracted number counts better represent 
the actual source population than the initial prior, making 
it possible to mitigate this bias through an iterative process 
that adjusts the prior based on the extracted counts. We ap- 
ply this iterative method to the AzTEC/COSMOS number 
counts estimate in the next section. 



2.3 AzTEC/COSMOS Number Counts 

The AzTEC/COSMOS integrated number counts are shown 
in Figure [3] This field shows an excess of sources at all 
fluxes when compared to the scaled SCUBA/SHADES re- 
sults (solid line). The number counts results are relatively 
insensitive to the choice of prior, with the initial analysis 
(filled circles) in agreement with those produced using an 
iterative prior (open squares) . The "robust" source criterion 
of P(Si < 0\S m ,&m) < 0.05 is somewhat more sensitive to 
the chosen prior, with the equivalent S/N threshold in a 
<r m = 1.3 mJy region being 4.2 and 4.0 for the initial and 
final iterative priors, respectively. The final iterative prior 
deems a larger number of sources as robust compared to the 
initial prior due to the number of sources lying in this S/N 
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Figure 3. AzTEC/COSMOS integral source counts derived from 
the most robust AzTEC sources in the field using the tech- 
niques described in £j |2,2l Filled circles (confidence bars omitted) 
represent the extracted counts using the initial assumed scaled 
SCUBA/SHADES prior (solid curve). Empty squares and 68% 
confidence bars represent the extracted counts after iteratively 
adjusting the prior to best represent the results of this survey. 
The dashed, dash-dotted, and dotted curves represent differential 
number counts fits to Eq. [2] with free parameters [S",7V*], [S'], 
and [TV*], respectively. The horizontal dashed line represents the 
"survey limit" , defined here as the source density that will Poisson 
deviate to zero sources (per 0.15 deg 2 ) 32.7% of the time. Open 
triangles and associated error bars repre sent the number c ounts 
found in the AzTEC GOODS-N survey llPerera et alj|2008h us- 
ing the same techniques and have been shifted 4% to the left for 
clarity. 
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range. The corresponding effect on the survey completeness 
keeps this from being a runaway process, with the iterative 
number counts quickly converging within a few iterations. 

The differential and integrated number counts of the 
AzTEC/COSMOS field are presented in Table □ We fit the 
differential number counts to Equation [2] using Levenberg- 
Marquardt minimisation, incorporating the data covariance 
matrix to account for correlations between flux bins. Vari- 
ous fits to the data are presented in Table [2] and are shown 
in Figure [3] Given the size and depth of this survey, we con- 
strain the parametric fits to flux bins between 2 and 10 mjy 
to avoid bins that are poorly sampled and prone to system- 
atic errors. This range of flux values is relatively insensi- 
tive to the power-law parameter a (Equation [2J; therefore, 
we fit the data while holding a constant at values of —2 
and — 1, which represent the SCUBA/SHADES result and a 
pure exponential, respectively. We also present similar fits to 
the SCUBA/SHADES number counts (jCoppin et alJl200d ) 
for comparison. The parametrised AzTEC/COSMOS re- 
sults provide the maximum constraint on differential source 
counts at fluxes ~ 4-5 mjy (depending on the parametrisa- 
tion). For example, a two parameter (S',N*; a fixed to —2) 
fit to Equation [2] constrains the AzTEC/COSMOS differ- 
ential counts at 4.5 mjy to 84 ± 17 deg~ 2 mjy -1 . 

Uncertainty in the flux calibration of the 
AzTEC/COSMOS survey is not included in these cal- 
culations and we believe it to be sub-dominant to the 
formal errors of the source flux, number counts, and 
fitted parameters. Calibration error estimates for individ- 
ual observations du ring this observing season are 6-13% 
( Wils on" et al.l [2008). Any normally distributed random 
component of this error will be reduced in the final co-added 
map since this survey is composed of multiple observations 
spanning many nights/weeks and calibration measure- 
ments. Systematic error in the calibration is believed to be 
dominated by the 5% uncer tainty in the flux dens ity of our 
primary calibrator, Uranus (|Griffln fc Orton]|l993h . 



Flux Density dN/dS Flux Density N(>S) 

(mjy) (mJy-Meg- 2 ) (mjy) (deg" 2 ) 



1.40 1 706 --1534. 1-0 2610.1^1^ 

2.42 535 -±292. 2 -° 904 -±3i9: 

3.43 204 --S2. 3 -° 369 -±l03. 

4.43 84 -±33. 4 -° 165 -±48. 

5.44 38. 5.0 81 --27. 
6.44 19-liJ: 6 -° 43 --it 
7.44 ll.+f; 7.0 24 -tll' 
8.44 6.5tgg 8.0 13 -±ib. 
9.44 3.9±3.6 9.0 & - 7 -6.7 



Table 1. AzTEC/COSMOS differential and integral number 
counts using iterative adjustment of the prior. The differential 
number counts flux bins are 1 mjy wide and span interger flux 
values, with effective bin centres (first column) weighted accord- 
ing to the assumed prior. The lowest flux bin listed is particularly 
sensitive to uncertainties in the prior and such systematics are not 
included in the given 68% confidence intervals. 



Data Set 


S' 


N* 


a 


x 2 




(mjy) 


(deg" 2 ) 






Az/COS 


1.83 ±0.41 


4420 ± 2720 


-2 


0.21 


Az/COS 


1.72 ±0.12 


5200 


-2 


0.28 


Az/COS 


1.36 


9610 ± 1970 


-2 


1.89 


SHADES 


3.36 ± 0.49 


5200 ± 1790 


-2 


0.23 


Az/COS 


1.31 ±0.23 


3570 ± 1790 


-1 


0.59 


SHADES 


2.39 ±0.27 


4370 ±1170 


-1 


0.21 


MODELx.lmm 


1.34 


5280 


-2 





Table 2. Fit results to the differential number counts and re- 
spective covariance matrix of th e AzTEC/COSMO S (1100 um) 
and SCUBA/SHADES (850 uim lCoppin et al.ll2006t ) surveys. All 
fits are to the modified Schechter function given in Equation [2] 
while holding various parameters constant (those with no uncer- 
tainty given). To avoid strong systematics at low flux, all fits are 
limited to data with Snoo^m ^ 2 mjy and Sssoum 4 mjy 
for the AzTEC and SCUBA surveys, respectively. The last row 
represents our assumed 1.1 mm blank-field model for the ini- 
tial prior (scaled SCUBA/SHADES), x 2 values are unrealistically 
low, likely due to a combination of: (a) our uncertainties being 
assumed as Gaussian in the fit; and (b) additional correlation 
not accounted for in the linear Pearson covariance matrix con- 
structed through the bootstrap sampling method (§ 12. 21 . These 
effects are also seen in the SCUBA/SHADES implementation of 
this algorithm jCoppin et alj|2006ft . 

2.4 Blank-Field Model 

It is immediately apparent that the AzTEC/COSMOS field 
is rich in bright sources when compared to other 1.1 mm 
surveys (see § I2.5|l . In order to quantify the significance 
of this potential over-density, we must first adopt an accu- 
rate characterisation of the true background ( "blank-field" ) 
population. The tightest published constraint on the SMC 
populat ion is provided by the 850 um SCUBA/SHADES 
survey l|Coppin et alJl200d ). which we convert to 1100 urn 
assuming an 850 (i.m/1100 urn power-law spectral index 
of 3.5. This scaling is roughly consistent with the inte- 
grated numb er counts of the 1. 1 mm Bolocam Lockman 
Hole survey ijLaurent et al.ll2005T h which partially overlaps 
with the SCUBA/SHADES survey. Assuming the SCUBA 
and AzTEC observations are in the Rayleigh- Jeans regime 
of optically-thin thermal dust emission from the SMGs, our 
scaling is consistent with the sub-mm spectral indexes of 
bright IR galaxies in the local universe (|Dunne et al ] |200d : 
iDunne fc EaleslbOOll ). 

Using a scaled version of the number counts measured 
at a different observation wavelength carries the inherent 
risk that the two bands are sensitive to significantly differ- 
ent (although overlapping) source populations, as evidenced 
by the possible ex istence of sub- millimetre drop-outs (SDOs, 
iGreve et al.ll2008h . The SCUBA 850 um surveys would be 
relatively insensitive to these proposed SDOs due to a com- 
bination of high redshift (z S> 3) and/or unusual spectral 
energy distributions (e.g. Td ust ~ 10 K). Therefore, it is 
important to verify the blank-field model with a direct mea- 
surement of 1.1 mm population. 

The most robust characterisation of the 
AzTEC/COSMOS over-density comes through com- 
parison to similar analyses of other AzTEC 1.1 mm surveys, 
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which eliminates systematics between different instruments 
and minimises those related to calibration. The best 
AzTEC 1.1 mm blank- field constraints are being provided 
by the AzTEC 0.5 deg 2 survey of the SHADES fields. 
Initial results of the AzTEC/SHADES survey (using nearly 
identical algorithms as those applied to AzTEC/COSMOS) 
are consistent with our scaling of the SCUBA/SHADES 
counts. Our number counts model falls in the higher regions 
of the AzTEC/SHADES uncertainty interval (modelled 
differential counts are roughly +0.5ct to +2. Oct above the 
average AzTEC/SHADES counts in the flux range explored 
here) and is within the field-to-field variations measured 
in those large surveys; therefore, we believe our model 
represents a conservatively high estimate of the blank-field 
counts that is appropriate for robust qualification of the 
potential over-density. 

We note that t he AzTEC survey of the GOODS-N field 
ijPerera et al. 1 120081 ) finds an SMC number density that is 
somewhat higher than our blank-field model; however, our 
model is within the ~ Ict uncertainty of that survey's inte- 
grated number counts (Figure [3] note that the data points 
are correlated). The AzTEC/GOODS-N results imply an 
S' parameter (S' — 1.25 ± 0.39 mjy) that is consistent 
with our general scaling of the SCUBA/SHADES counts, 
but suggest systematically higher number counts (i.e. larger 
TV* parameter). The small size of the GOODS-N survey 
(0.068 deg 2 ) makes it highly susceptible to cosmic variance 
and clustering, thus reducing it's viability as a measurement 
of the average sky. It also does not significantly constrain 
the bright (S > 5 mjy) 1.1 mm source counts where the 
AzTEC/COSMOS over-density is most apparent (§ [23)l . 
The SCUBA survey of GOODS-N (|Borvs et al] 12003 ) al- 
ready suggests the field may be overly rich in sub-millimetre 
sources, with number counts systemati cally higher than see n 
in the SCUBA/SHADES blank-field (|Coppin et al.l l2006h . 
although the analyses of these two SCUBA surveys differ 
significantly and the difference in number counts could be 
partially systematic. 

2.5 SMG Over-density 

The source catalogue presented in Paper I suggests that the 
AzTEC/COSMOS field has a significantly larger number of 
bright 1.1 mm sources than might otherwise be expected for 
a survey of this size and depth. The density of sources in the 
AzTEC/COSMOS field with raw measured fluxes > 6 mjy 
is 3 times higher (14 sources in 0.15 sq. deg. field) than in 
the 1.1 mm Bolocam Lock man Hole survey of s imilar depth 
(3 sources in 0.09 sq. deg.; lLaurent et aill2005l ). The seven 
brightest AzTEC sources in the COSMOS field have been 
imaged using the Submillimeter Array (SMA) at 890 um, 
and they are sho wn to be single, un resolved sources at 2 
arcsec resolution (|Younger et al.ll2007T ). 

We compare the AzTEC/COSMOS number counts to 
the blank-field model discussed in § 12.41 Figure [3] shows that 
the AzTEC/COSMOS integrated source count estimates are 
clearly in excess of the scaled SCUBA/SHADES counts. 

To estimate the probability of this excess happening 
by chance, we compare the number of robust sources de- 
tected in the AzTEC/COSMOS survey to the number re- 
covered in simulated maps. In Figure [4] we show the distri- 
bution of the number of recovered sources in 10,000 simula- 



tions, each populated with a random realisation of the scaled 
SCUBA/SHADES counts. On average, 12.1 sources are re- 
covered from each of the simulated maps, which is in agree- 
ment with the semi-analytical expectation value of 11.2 (cal- 
culated from the scaled SCUBA/SHADES results and simu- 
lated completeness of this survey) and the expected number 
of false detections ({NfaUe} ~ 1.2). Application of the same 
source criteria (5% null threshold, scaled SCUBA/SHADES 
prior) to the real map results in 23 robust sources (32 if 
using the iterative prior), which is greater than in 99.7% of 
the simulations. The AzTEC/COSMOS source over-density 
is even more significant in the number of very bright sources, 
with 11 detected at S/N ^ 5 (S m £ 6.2 mjy). Ten- 
thousand simulations of the blank-field model could produce 
no more than 6 such detections in a single map, thus infer- 
ring a S> 4ct significance in the number of bright sources. 

Assuming the scaling of the SCUBA/SHADES counts 
accurately represents the blank-field SMG population at 
1.1 mm, the parametric fits shown in Table [5] favour the 
over-density being described as a shift in the flux pa- 
rameter S" over an increase in the normalisation param- 
eter N* . This is consistent with an apparent over-density 
caused by uniform amplification of the source fluxes. How- 
ever, this solution is degenerate with an alternative scal- 
ing of the SCUBA/SHADES results - scaling with a flat- 
ter spectral index of ~ 2.6 also produces a good fit to 
the AzTEC/COSMOS number counts (dash-dotted curve 
of Figure [3]). Therefore, an accurate representation of the 
blank-field population (§ I2.4|l is critical to properly quanti- 
fying the over-density. 

If taken alone, the over-density of SMGs in the 




Number of Sources 

Figure 4. Distribution of the number of robust sources detected 
in 10,000 simulations (using 300 unique noise realisations) of 
the AzTEC/COSMOS observations when randomly populating 
the astronomical sky with the scaled SCUBA/SHADES number 
counts (solid histogram). Using the same source criteria, there 
are 23 robust sources detected in the AzTEC/COSMOS map 
(dashed vertical), which is greater than 99.7% of the simulations. 
A Gaussian fit to the simulation results (thin solid curve) shows 
23 sources to be a 3.3o" outlier. The difference between the sim- 
ulation mean (solid vertical) and the semi-analytic expectation 
value (dotted vertical) reflects the number of false detections (i.e. 
~ 1 per map) for the scaled SCUBA/SHADES prior and chosen 
source threshold. 
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AzTEC/COSMOS field would likely be explained away as 
simple cosmic variance in the SMG population as traced 
in a 0.15 deg 2 field. However, in the following sections we 
demonstrate that the over-density is due, in part, to fore- 
ground structure in the COSMOS field. Only with the rich 
multi- wavelength coverage of the COSMOS field and the rel- 
atively large size of the AzTEC map is this analysis possible. 



3 CORRELATION BETWEEN AZTEC 
SOURCES AND LARGE SCALE 
STRUCTURE IN THE COSMOS FIELD 

Having shown that the AzTEC/COSMOS field exhibits 
a significant excess of bright SMGs with respect to our 
adopted blank-field model, we explore the possibility that 
this over-density is due, in part, to a correlation of 
AzTEC sources with the prominent large-scale structures 
at z ;$ 1.1 identifie d in this portion of the COSMOS field 
|Scoville et al.l l2007s) . All correlation tests in this section are 
limited to the inner 0.15 deg 2 region of the AzTEC map 
where the uniformity in coverage simplifies the analysis. 

The smoothed gal axy density map pr oduced by the 
COSMOS consortium (IScoville et al.1 120071 . see Figure [5} 
shows a collection of dense regions in the the AzTEC cov- 
ered area. We first look for coincidence with AzTEC sources 
by cross-correlating the surface-density of optical-IR galax- 
ies in this map with the AzTEC source positions using a 
bi-d imensional Ko l mogorov-Smirnov (K-S) test of similar- 
ity (|Peacockl 1 19831 : iFasano fc Franceschinil Il987l ). Restrict- 
ing this analysis to the 30 robust AzTEC sources detected 
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Figure 5. Smoothed surface density map of ga laxies derived 
from the optical-IR catalogue of COSMOS galaxies fSc oville et al.l 
l2007t) in the 0.15 sq. deg area surveyed by AzTEC, where darker 
colours indicate more densely populated areas of the sky. The 
cross and plus symbols represent AzTEC sources detected at 
signa l-to-noise ratios S/N ^ 4 and 4 > S/N > 3.5 ( Sco tt et all 
2008), respectively. 



at a S/N ^ 4, wh ich have an estim ated false detection rate 
of less than 7% l|Scott et al.ll2008ft . we find that the test 
cannot reject the null hypothesis that the distribution of 
AzTEC-source positions follows the surface density distri- 
bution of optical-IR galaxies. The test concludes that the 
difference between the SMG and optical-IR populations is 
smaller than 93.7% of the differences expected at random 
due to sample variance, often referred to as rejecting the 
null hypothesis at the 6.3% level, thus suggesting the distri- 
butions could indeed be similar. 

The significance of the SMG positional correlation with 
the large-scale structure is further quantified by compar- 
ing the K-S D-statistic of the SMG catalogue to that of 
a homogeneous random distribution of the same number 
of sources, under the null hypothesis that they follow the 
surface density of optical-IR galaxies. This test determines 
that the AzTEC/COSMOS source distribution follows the 
optical-IR distribution more strongly than 98.9% (~ 2.5a) 
of the random-position catalogues. The result is somewhat 
less significant, 91.1%, if we expand the comparison to the 
full S/N > 3.5 AzTEC/COSMOS catalogue, which is likely 
due to the increased number of false detections (from ~ 1 
to ~ 11) at this lower S/N threshold. These false detections 
(noise peaks) are inherently random and homogeneous in 
their distribution and dilute the correlation signal. 

It is possible that only a fraction of the AzTEC source 
positions are correlated with the prominent large-scale struc- 
tures detected in the COSMOS galaxy density map while 
a subset of randomly distributed source positions dilutes 
the sensitivity of the quadrant-based bi-dimensional K-S 
statistic discussed above. Therefore, we further test the hy- 
pothesised correlation by comparing the surface density of 
optical-IR galaxies within a small area surrounding AzTEC 
positions to that surrounding random positions in the map. 
Figure [6] shows the distribution of the galaxy densities at 
z 1.1 projected within 30" (1.7 pixels in the smooth 
galaxy density map) of the AzTEC source positions, com- 
pared to the galaxy densities found around random positions 
within the AzTEC survey region. The two distributions are 
clearly different, with a one-dimensional K-S test rejecting 
the null hypothesis of identity at > 99.99% and 97.2% lev- 
els for the S/N 4 and S/N ^ 3.5 catalogues, respec- 
tively. The mean number of nearby optical-IR galaxies at 
S/N > 4 (S/N > 3.5) AzTEC source positions is larger 
than that at random positions in the map at a significance 
of 99.99% (99.5%) according to the non-parametric Mann- 
Whitney (MW) [/-test. 

We can search in redshift space for the structures that 
contribute the most to the coincidence between AzTEC 
sources and the galaxy density in their "line of sight" 
using the photometric redshifts of the optical-IR popula- 
tion (|llbert et al.l 120081 ) . which have a mean accuracy of 
\Az\/(l+z) » 0.01-0.02. Figure[7]shows a bar-representation 
of the MW probabilities that the mean integrated galaxy 
density around AzTEC sources is significantly larger than 
that around random positions in the map for various red- 
shift slices. There is positive signal (J> 2a) arising at differ- 
ent redshift slices, most notably at z ~ 0.65. At redshifts 
2 > 1.1, the number of galaxies detected at optical-IR wave- 
lengths decreases significantly, and the level of correlation 
found with AzTEC sources is well below the 2<r threshold. 

The most prominent contribution to the AzTEC- 
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normalized galaxy dersity 



Figure 6. Histogram of the galaxy density at z <J 1.1 found 
within 30 arcsec of (a) AzTEC/COSMOS source candidates de- 
tected at S/N ^ 3.5 (dashed line) and (b) random positions in 
the AzTEC mapped area of COSMOS (solid line). The two pop- 
ulations are different at the 97.2% confidence level, using a one- 
dimensional K-S test. The galaxy densities are normalised to the 
mean galaxy density in the full AzTEC-covered area. 




Figure 7. Bar-representation of the Mann- Whitney probability 
that the mean galaxy density around AzTEC sources at a given 
redshift slice is significantly larger than the mean galaxy den- 
sity around random positions. Horizontal dotted lines represent 
the la, 2a, and 3a significance levels, respectively. The blue- 
dotted curve shows the relative number of optical/IR galaxies 
contained within each redshift slice within the AzTEC covered 
area (AT gal ( 2 )/AT total * 10. + 0.60). 



optical/IR correlation lies at 0.6 z ^ 0.67, with the red- 
shift slices within this range having MW probabilities of 
difference up to 99.98%. The smoothed galaxy density map 
for this redshift range is shown in Figure [8] Two promi- 
nent large-scale s tructures have been identified (Structures 
#1 and #24 in IScoville et ail 120071 ) within this redshift 
slice. Structure #1 at z = 0.73 ± 0.27 has 1767 optical- 
IR galaxy members and approximately spans (FWHM) 
A RA = 0.22 deg and A Dec = 0.17 deg. Structure #24, 
a less massive but very compact system that is X-ray de- 
tected, has 85 galaxy members and is at z ~ 0.61. Structure 
#24, however, does not appear to contribute to the correla- 
tion, as no AzTEC sources fall within its primary extension. 
Structure #1 contains a rich core and represents a massive 



cluster (~ 10 Mq) at z ~ 0.73, which is c learly seen in 
the C OSMOS weak-lensing c onvergence map |Massev et al.l 
120071 ) and in X-ray emission (|Guzzo et al.ll2007l ). This clus- 
ter lies outside the redshift span of strong correlation, but 
the filamentary structure that leads to it is part of the red- 
shift slice under analysis (see Figure |SJ). 

We next assess whether the substructures within Struc- 
ture #1 are the main contributors to the observed corre- 
lation. Figure [9] shows the distribution of galaxy densities 
around AzTEC sources and around random positions in the 
collapsed 0.6 <, z <J 0.67 map. The means differ at the 
99.8% confidence level according to the MW (7-test. If we ex- 
clude a circular region around the cluster centre with radius 
9 = 1.5'(~ 0.6 Mpc), which contains both the cluster-core 
and the clu ster-outskirt region s seen by the X-ray tempera- 
ture profile l|Guzzo et al.ll200"7l ). the significance of the differ- 
ence is 99.94%. Excluding a larger circular region of radius 
8 = 6'(~ 2.4 Mpc), which represents the FWHM of the full 
Structure #1, the MW significance decreases to only 98.6%. 
This demonstrates that although AzTEC sources do corre- 
late with the galaxy densities associated with the extended 
Structure #1, the less prominent large-scale structure across 
the rest of the map is also well-correlated with the AzTEC 
positions. 

Figure [lOl shows the optical/IR galaxy density map for 
the redshift slice 0.24 <, z <J 0.26, which is also a large con- 
tributor to the overall correlation between the large-scale- 
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Figure 8. Smoothed surface density map of galaxies at 0.60 <J 
z <: 0.6 7 detected at optic al-IR wavelengths by the COSMOS 
survey llScoville et al.l l2007t) ■ The large-sca le structure at z = 
0.73 ± 0.26 detected bv lScoville et al.l fcoOTl) is marked as Struc- 
ture #1 and the large circle (6' diameter). This large-scale struc- 
ture has a peak over-density at z ~ 0.73, outside of the red- 
shift range of this figure, and is identified as a massive cluster 
jGuzzo et al.ll2007fl . The yellow circle (1.5' diameter) marks the 
spatial extent of this cluster as traced by the X-ray contours. 
Another rich cluster, at z ~ 0.61, is marked as Structure #24. 
Symbols are the same as in Figure [5] 
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Figure 9. Histogram of the fraction of optical-IR selected galax- 
ies at 0.60 z < 0.67 found within 30 arcsec of AzTEC source 
candidates (red dashed line), and around random positions within 
the AzTEC-mapped area of COSMOS (black solid line). The 
blue dotted-line histogram represents the number of optical-IR 
galaxies around AzTEC sources, excluding the single AzTEC 
source that falls within the X-ray trac ed cluster-outskirt region 
(0 <y 1.5 arcmin from the cluster centre. fCuzzo et al. 2007), while 
the green dash-dotted histogram excludes the full 6 arcmin radial 
structure identified as Structure #1. The mean values of these his- 
tograms are represented at the top of the figure as vertical bars. 
The galaxy densities are normalised to the mean galaxy density 
in the full AzTEC-covered area such that the mean density of 
random positions is 1. 



structure of the field and AzT E C sou rces (Figure [TJ. Struc- 
ture #22 from IScoville et ail l|2007t ). with ~ 67 possible 
galaxy members at z ~ 0.26±0.11, is the main cluster in this 
redshift slice and is also detected in X-ray. However, as with 
the portion of Structure #1 in the 0.60 < z < 0.67 slice, 
this system does not dominate the overall correlation with 
AzTEC sources: the mean galaxy density around AzTEC 
sources differs from random locations at the 99.1% level af- 
ter exclusion of the A RA ~ 0.06 deg and A Dec ~ 0.14 deg 
area of influence of the cluster. Similarly, we find that the 
prominent contributions of other redshift slices (e.g. z ~ 0.33 
and z ~ 0.8) to the overall spatial correlation are not due to 
single compact structures. 

It appears that the observed correlations are not dom- 
inated by the clusters in the field, thus it is not surprising 
that the AzTEC positions are, in general, l ess correlated 
with the weak-lensing mass map of COSMOS l|Massev et al.l 
2007), which is particularly sensitive to the most massive 
structures like the z ~ 0.73 cluster (see Figure ITT)) . The null 
hypothesis that the distribution of masses found within 30" 
of AzTEC positions is the same as that found around ran- 
dom positions in the weak-lensing map is "rejected" at only 
the 60% level (K-S test), and their means differ at the 91.5% 
level (MW [/-test). 



4 DISCUSSION 

The distribution of AzTEC sources is correlated with the 
large-scale distribution of optical/IR galaxies at z <J 1.1 and 
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Figure 10. Smoothed surface density map of galaxies at 0.24 ;C 
z <J 0.26 detected at optical-IR wavelengt hs by the C OSMOS sur- 
vey, which includes Structure #22 of lScoville et all ||2007|) . Sym- 
bols are the same as in Figure [5] 



the primary (but not unique) contributors to this signal are 
located at redshifts 0.60 <; z & 0.65 and 0.24 <, z <, 0.26 
(section §03). The correlations in these redshift regimes are 
robust, with mean optical/IR densities differing from ran- 
dom distributions at the 99.1 to 99.98% level (MW [/-test). 
For the seven AzTEC so urces that have been followed up 
with SMA interferometry l|Younger et al. l2007h and the ad - 



ditional 14 that have radio detections ( Scott et al.ll2008h . 
secure optical/IR counterparts have been identified. The 
optical-IR and FIR-mm-radio photometric redshifts of these 
sources place the majority of these objects at z J> 3; there- 
fore, AzTEC detected sources are most likely background 
systems to the z <, 1.1 galaxy densities shown in Figure [5] 
This is not surprising, given that the population of SCUBA 
SM Gs with radio counterparts have a median redshift of 
2.2 dChapman et alJliool . 120051 : lAretxaea et"al]|2003l . 120071 : 
iPope et al.ll2005l ). 

The amplification caused by massive clusters at inter- 
mediate redshifts (z ~ 0.2 — 0.4) has been used to detect and 
study the sub-milli metre galaxy p opulation since the first 
SMG surveys (e.g. ISmail et al.lll997l ). Lensing is expected 
to occur also in and around the z w 0.73 cluster detected in 
the COSMOS field, but only 4 of the 50 S/N > 3.5 AzTEC 
source-candidates are projected within 2' of the dense clus- 
ter core. Removing these sources/regions from the analysis 
of § [3] has little effect on the correlation strength. Thus the 
correlations seem to be tied to the general z ^ 1.1 large-scale 
structure in the field. The same result holds if we exclude 
the other prominent structures in this field, Structures #24 
at z « 0.61 and #22 at z ~ 0.26. Furthermore, the bulk of 
AzTEC sources do not significantly correlate with the weak- 
lensing map, which is particularly sensitive to the mass con- 
tained in rich clusters. 
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Lensing of the sub-mm galaxy population by foreground 
low-redshift structures has been claimed in the correla- 
tion analysis of 39 sub-millimetre galaxies detected in 3 
disjoint fields with the density of R < 23 mag galaxies 
(Al maini et al.ll2005h . which statistically lie at (z) ~ 0.5. 
It was argued that the bright Ssso/jm > 10 mjy sources are 
found to c luster preferen ti ally a round the highest-density 
areas, and lAlmaini et al.l (|2005h estimate that 20-25% of 
the sub-millimetre galaxy population is subject to lensing 
by foreground structures. We no te that a similar s tudy per- 
formed in the GOODS-N region |Blake et alj2006j ) found no 
detectable correlation between 35 SCUBA-selected SMGs 
and the optically-selected galaxy populations at z ^ 0.8. 
This difference in correlation strength may be related to cos- 
mic variance of foreground structure on the scale of these 
maps. There also exist potential cases of lensing by in- 
dividual galaxies, with some SCUBA sources being incor- 
rectly identified as low-redshift galaxies due to intervening 
foreground galaxies that lie directly alon g the line of sight 
[|Chapman et al.ll2002l ; iDunlop et afll2004h . Since it includes 
a high-density region within the COSMOS field, the AzTEC 
survey is sensitive to all of these types of amplification, and 
we have demonstrated that there is a positive correlation 
with the large-scale structure. Inspection of the optical/IR 
counterparts of the 21 AzTEC galaxies with radio and/or 
sub-mm interferometric positional accuracy, including the 
7 sources known to have sub-m illimetre emission on scales 
6 < 1.2" (| Younger et al.ll2007h . show no obvious signs of 
strong galaxy-galaxy lensing and hence any amplification of 
this sub-sample must be attributed to weak-lensing. 

If our blank-field number counts model (§ I2.4|l accu- 
rately represents the intrinsic (non-amplified) SMC popula- 
tion in the AzTEC/COSMOS field, then the observed num- 
ber density of sources is also consistent with weak-lensing of 
the background SMC population. Parametric fits to the flux- 
corrected number counts (§ 12.31 see also Table [2] and Fig- 
ure [3]) show that the relative over-density of sources can be 
fully explained as a systematic increase in the parameter S' , 
which is consistent with an average flux amplification (e.g. 
lensing) of the source population by ~ 30%. Conversely, the 
number counts data are only marginally consistent with a 
simple increase in the normalisation parameter N* , thus dis- 
favouring a uniform physical over-density (e.g. cosmic vari- 
ance) of sources in this field as the sole cause of the observed 
over-density. Additionally, any over-density due to variance 
and/or clustering can not explain the correlation of AzTEC 
sources to the z ^ 1.1 structure, as the AzTEC sources are 
likely background sources and not physically associated with 
the z ^ 1.1 structure. 

An alternative cause of the number counts over-density 
can be imagined as an additive flux source (e.g. dense screen 
of faint foreground sources) confused with the blank-field 
sources. However, the AzTEC/COSMOS map has been fil- 
tered for point s ource detection and has a mean of zero 
l)Scott et al.ll2008T ) , which leaves the map insensitive to high- 
density or uniform millimetre flux sources that span large 
spatial scales. Furthermore, the positions of AzTEC sources 
are not strongly correlated with the most dense and com- 
pact foreground regions (i.e. clusters) that could otherwise 
be potential sources of additional mm-wave flux in our map 
(e.g. the Sunyaev-Zel'dovich effect). 

The significance of the spatial correlation between 



AzTEC sources and the intervening large-scale structure 
contrasts with the lack of a similar detectable signal a mong 
the sources discovered by COSBO (jBertoldi et al.ll2007l ). the 
1.2 mm MAMBO survey to the south and adjacent to the 
AzTEC surveyed area (see Figure [TTj) . If we repeat the anal- 
ysis performed in §|3]with the MAMBO catalogue, we do not 
find a significant correlation with the COSMOS optical/IR 
galaxy surface density; the probability that the galaxy densi- 
ties around MAMBO sources are different from that around 
random positions in the map is only 87% (K-S test). This 
lack of a significant correlation signal may be due in part to 
the smaller catalogue of significant sources in the COSBO 
field and the overall lack of significant foreground structure 
in much of the COSBO covered area. 

The association between COSBO sources and the weak- 
lensing derived mass-map, however, is stronger with a 99.4% 
probability (K-S test) that the distribution of mass around 
MAMBO source locations is different from that of random 
positions in the COSBO survey region. This signal is domi- 
nated by a group of ~ 7 — 8 of the most significant MAMBO 
sources close to two compact mass spikes, which are identi- 
fied with X-ray bright over-densities consistin g of a total of 
~ 12 7 galaxies at z « 0.24 (Structure #17 in Scoville et al.l 
and are likely clusters. The possible spatial correla- 
tion of COSBO sources with foreground structures, there- 
fore, might be of a somewhat different nature than that of 
AzTEC sources. The COSBO region may be witnessing am- 
plification caused by the two clusters revealed by the weak- 
lensing map, while the AzTEC sources are more likely am- 
plified by galaxies contained within the more tenuous fil- 
amentary large-scale structure, which is so tenuous in the 
COSBO field that it provides no significant signal in the 
sample of MAMBO sources. 



5 CONCLUSIONS 

The central 0.15 deg 2 of the AzTEC/COSMOS survey 
shows a significant over-density of bright 1.1 mm-detected 
SMGs when compared to the background population in- 
ferred by other surveys. We find that this over-density 
cannot be explained as sample variance of the blank-field 
SMG population. The SMG positions are significantly cor- 
related with the z 1.1 optical/IR galaxy density on the 
sky, which is believed to be in the foreground of nearly all 
AzTEC/COSMOS SMGs. Both the spatial correlation and 
the AzTEC/COSMOS SMG number counts are consistent 
with gravitational amplification of the blank-field SMG pop- 
ulation. The lack of st r ong co rrelation to the weak-lensing 
maps of iMassev et al.l l|2007h indicates that this amplifi- 
cation is primarily due to weak-lensing by the large-scale 
structure as opposed to lensing by the compact and massive 
clusters in the field. SMGs detected in a different part of 
the COSMOS field by the 1.2 mm COSBO survey are also 
spatially correlated to the z 1.1 structure, however, this 
correlation is dominated by two compact structures (likely 
clusters) in the field. The lack of significant large-scale struc- 
ture (i.e. lensing opportunities) in the rest of the COSBO 
survey region results in COS BO number counts t hat are con- 
sistent with the blank-field l|Bertoldi et al.ll2007l ) - a strong 
contrast to the significant SMG over-density and rich fore- 
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Figure 11. Weak-lensing convergence mass map jMassev et al.l 
120071) of the AzTEC and MAMBO surveyed regions of COSMOS. 
Crosses mark the millimetre source positions from the two cata- 
logues and follow the notation of Figure [5] 



ground structure found in the nearby AzTEC/COSMOS 
field. 
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